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Abstract. Aerosol concentrations and size distributions in the middle and upper troposphere 
over the remote Pacific Ocean were measured with a forward scattering spectrometer probe 
(FSSP) on the NASA DC-8 aircraft during NASA’s Global Backscatter Experiment (GLOBE) in 
May-June 1990. The FSSP size channels were recalibrated based on refractive index estimates 
from flight-level aerosol volatility measurements with a collocated laser optical particle counter 
(LOPC). The recalibrated FSSP size distributions were averaged over 100-s intervals, fitted with 
lognormal distributions and used to calculate aerosol backscatter coefficients at selected wave- 
lengths. The FSSP-derived backscatter estimates were averaged over 300-s intervals to reduce 
large random fluctuations. The smoothed FSSP aerosol backscatter coefficients were then com- 
pared with LOPC-derived backscatter values and with backscatter measured at or near flight level 
from four lidar systems operating at 0.53, 1.06, 9.1 1, 9.25, and 10.59 /xm. Agreement between 
FSSP-derived and lidar- measured backscatter was generally best at flight level in homogeneous 
aerosol fields and at high backscatter values. FSSP data often underestimated low backscatter 
values especially at the longer wavelengths due to poor counting statistics for larger particles 
(>0.8 jum diameter) that usually dominate aerosol backscatter at these wavelengths. FSSP data 
also underestimated backscatter at shorter wavelengths when particles smaller than the FSSP 
lower cutoff diameter (0.35 jx m) made significant contributions to the total backscatter. 


1. Introduction 

Tropospheric aerosols which lie above the planetary 
boundary layer (PBL) have in the past 15 years become an im- 
portant subject of research for space-borne measurements of 
winds in this region of the atmosphere. Currently only sparse 
data are available on tropospheric winds over remote ocean 
areas and passively borne aerosols are a surrogate to determine 
wind fields over this region using a space-based Doppler lidar 
system [ Menzies , 1986; Baker et ai, 1995]. However, the 
aerosol properties useful for wind measurements in these 
regions have not been extensively studied. Rothermel et ai 
[1989; 1995a] reported a low background backscatter mode in 
the middle and upper troposphere based on statistical reanaly- 
sis of infrared ground-based lidar climatological measurements 


institute for Global Change Research and Education and Earth 
System Science Laboratory, University of Alabama, Huntsville, 
Alabama. 

2 NASA Ames Research Center, Moffett Field, California. 

^Institute for Global Change Research and Education and 
Universities Space Research Association, Huntsville, Alabama. 

4 School of Oceanography, University of Hawaii, Honolulu, Hawaii. 
■’NASA Marshall Space Flight Center, Huntsville. Alabama. 

6 NASA Goddard Space Flight Center, Greenbelt, Maryland. 

7 Jet Propulsion Laboratory, Pasadena, California. 

Copyright 1996 by the American Geophysical Union. 

Paper number 95JD03423. 

0 1 48 -0227/96/95 JD-03423$05. 00 


made over a number of years over Colorado and airborne 
Doppler lidar measurements over the United Kingdom. The 
weak magnitude of this background aerosol mode would pose a 
stringent test on the effectiveness of any satellite-borne 
Doppler laser wind sensor system to provide tropospheric 
wind data (D. A. Bowdle et al., unpublished data, 1995). 

Instrumentation which measures aerosol properties is 
diverse. Several types of sensors are required to measure unam- 
biguously a synergistic set of aerosol properties [ Prospero et 
al . , 1983]. For example, pulsed lidar systems can measure 
aerosol backscatter coefficient above or below flight level but 
only provide limited data on the aerosol size distributions if 
several wavelengths are used simultaneously. Detailed, real- 
time aerosol size distribution data are obtained primarily using 
optical particle counters, which can be used with measured 
aerosol composition to predict aerosol backscatter at flight 
level. When different sensors are used in this way, direct com- 
parison with measured aerosol backscatter data is required for 
corroboration, particularly where low aerosol concentrations 
prevail. To obtain data to guide design and performance studies 
of prospective satelite Doppler wind lidar, NASA conducted 
the Global Backscatter Experiment (GLOBE) survey mission 
flights in November-December 1989 (GLOBE I) and May-June 
1990 (GLOBE II) to sample the free troposphere aerosol in the 
Pacific basin [Bowdle et al. y 1991]. This area was selected 
because of the possible extensive occurrence of the back- 
ground aerosol backscatter mode in the midtropospheric 
region. The GLOBE II instrumentation aboard the NASA DC-8 
aircraft consisted of optical particle counters sampling at 
flight level and several lidar systems to measure aerosol 
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GLOBE 2 FLIGHT TRACKS, 1990 



Figure 1. Flight paths for all transit flights made during 
NASA’s GLOBE II May-June 1990 Pacific survey mission. 

backscatter coefficient at or near flight level for several 
wavelengths. 

This paper will report on selected aspects of the aerosol 
backscatter coefficient (ft, m' 1 sr* J ) data obtained in the second 
survey mission of GLOBE. In particular, the comparison is 
reported between backscatter coefficient data computed from 
NASA Ames Research Center (ARC) forward scattering spec- 
trometer probe (FSSP) size distribution data and that measured 
with four lidar systems operating at 0.53-, 1.06-, 9.1 1-, 
9.25-, and 10.59-/im wavelengths. Discrepancies are ex- 
plained by instrument performance characteristics and by spa- 
tial variability in aerosol microphysics and composition at or 
near Bight level. The results are used to determine the finest 


temporal resolution that can be achieved in a wide variety of 
atmospheric conditions for aerosol backscatter values derived 
from FSSP measurements. 

2. Experimental Details 

All data reported here were taken aboard the NASA DC-8 air- 
craft on the GLOBE II Pacific basin survey mission during 
May-June 1990 [Bowdle et aL , 1991]. Figure 1 shows the 
transit Bight tracks Oown, most of which were made near alti- 
tudes of 8 and 12 km. Several local Bights also were made and 
these generally consisted of stepped ascents and descents. In 
this paper the analysis will concentrate on (1) data from two 
case studies (transit Bights 6 and 14) and (2) combined data 
from Bights 6, 8, 10, 12, and 14. GLOBE instruments can be 
grouped into two classes: those measuring aerosol backscatter 
coefficient (lidars) and those measuring selected aerosol 
microphysical properties (optical particle counters) is the 
volume angular scattering at 180° from a particulate medium 
per unit length of the irradiating beam, per steradian centered 
about the direction of the incident beam). These instruments 
are brieBy described in the following sections. Table 1 shows 
the typical sample volume for each sensor described when the 
aircraft is Bying at 200 m s' 1 with an integration time of 
100 s and its relationship to the aircraft. 

2.1. Optical Particle Counters 

Optical particle counters were the primary instruments used 
on the DC-8 aircraft to measure the concentrations and size 
distributions of aerosols at Bight level. The FSSP-300 probe, 
which provided data for the main analysis in this paper, was 
manufactured by Particle Measuring Systems, Incorporated 
(PMS). The probe is designed to measure in situ size distribu- 
tions with fine temporal resolution in linearly distributed size 
bins over a 0.35- to 20 - fim diameter range [ Baumgardner et aL , 
1992; Pueschei et aL , 1994]. The FSSP pod unit was mounted 
on the extremity of the right wing. In-Bight data (counts/bin) 
were recorded at 1-s intervals and later converted to size- 
distributed concentrations averaged over 10-s intervals. The 
instrument was calibrated for size bin accuracy using 
monodisperse latex spheres. Sampling efficiency was assumed 
to be unity. Baumgardner et aL [1992] have estimated the 


Table 1. Comparison of GLOBE Sensor Sample Volumes 


Sensor 

— 

Laser 

Wave- 

length, 

/xm 

Lidar 

Pulse 

Rate, 

Hz 

Lidar Vertical 
Integration, 
m 

Beam 

Area, 

m 2 

Sample Volume 
Relative to A/c 

Sample 

Volume, 

m 3 

JPL lidar 

9.25 

4 

150 

1 

nadir or zenith 

6 x I0 4 

GSFC lidar 

0.53 

50 

75 

i 

nadir or zenith 

4 x 10 5 

MSFC CW 

1.06 
1 .54 
9.1 1 



3 x 10' 2 * 

50 m from RHS 

6 x I0 2 

lidars 1 

FSSP f 

10.59 

. 

- 

5 x lO 8 

of A/c 

RHS wing pod 

1 x 10- 3 

LOPC^ 

- 

- 

- 

1 x 10- 6 

A/c entry nozzle 

1.7 x lO’ 3 


GLOBE, 100-s integration time is used for all sensors. 

* Assumes a 5-m lidar beam focus depth, 50 m from the aircraft. 
'"Calculations based on an aircraft speed of 200 ms 1 . 

+A Bow rate of 17 cnT s’ 1 with virtual impactor. 
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Table 2. Threshold Levels and Uncertainties for Lidar Measured /? 


Lidar 

Wavelength, jum 

Threshold ft Level, 
m' 1 sr 1 

P Uncertainty 

GSFC pulsed 

0.53 

5xl0- 9 

±lx I0' 8 

GSFC pulsed 

1.06 

IxlO' 9 

±1x1 0' 9 

MSFC CW 

9.1 1 

6x1 O' 12 

±0.72xl0' 12 ’ 

JPL pulsed 

9.25 

IxlO' 1 1 at 2 km 
from A/c 

±2x10'' 1 

MSFC CW 

10.59 

2x10" 

±0.24x1 0' 1 '* 


*Using surface acoustic wave (SAW) signal processing. 


uncertainty in the total number concentration measured by the 
FSSP as 25% assuming a sampling error of 10%. 

The second system, a thermally preconditioned laser optical 
particle counter (LOPC) developed by the University of Hawaii 
[ Clarke , 1991] provided size-segregrated information on the 
aerosol physicochemistry in near real time. A PMS laser 
aerosol spectrometer (model LAS-X) was modified to provide 
256 channels of equal logarithmic width for particle diameters 
from 0.15 to 7.0 fi m. A virtual impactor (VI) upstream of the 
LOPC is employed to increase the concentrations of the larger 
particles greater than 0.5 fim diameter. The preheater system 


interfaced with the LOPC, in turn, heated the aerosol to 40°C, 
150 C, and 300 C. This allowed size-resolved measurements 
of sulfuric acid, ammonium sulfate, ammonium bisulfate, and 
refractory (e.g., dust, soot) aerosols to be obtained based on 
their volatility at these temperatures. Aerosol for sampling by 
the LOPC was bought into the aircraft through a blunt entry 
nozzle diffuser. Flow rate was maintained at isokinetic values 
based on measured flows and corrected for tube pressure, tem- 
perature, and large-particle loss [Porter et al . , 1992]. LOPC 
data were obtained in flight by alternately directing the sample 
flow through each preheater channel without the VI and 
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Plate 1 . False color cross section of \og(dnld log (D p )) for Bight 14 derived from coregistered (a) 20-s 
forward scattering spectrometer probe (FSSP) recalibrated and (b) laser optical particle counter (with virtual 
impactor) sensor data. Line at log )0 (0.39) in Plate lb corresponds to lower diameter limit for the recalibrated 
FSSP sensor. Note that diameter range (0.15 to 20 //m) and dnld log (D p ) range (10* 6 to 10 4 cm' 3 ) are both in 
logarithmic units. In Plate la the bar at the bottom gives the integrated particle numbers over 0.35-20 /im, 
while the two bars in Plate lb cover size ranges 0.16-7/im and 0.39-7 /im, respectively. 
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through the 40° C and 300°C channels with the VI. Several 
complete cycles were performed during a typical 10- to 15-min 
integration period to avoid confusing temporal variability 
with aerosol volatility. These sample periods were usually 
divided equally among the five preconditioning channels, so 
the total sample time in each channel was about 2-3 min. The 
LOPC was calibrated for sizing accuracy and constituent dis- 
crimination using well-characterized, laboratory-generated 
monodisperse particles with compositions similar to atmo- 
spheric aerosols. 

2.2. Lidars 

The lidar instrument complement to measure aerosol 
backscatter coefficients directly consisted of four lidar sys- 
tems operating at six wavelengths. Each of these systems (two 
pulsed and two continuous) will be briefly described. Table 2 
summarizes the threshold level and uncertainty in £ measured 
by each lidar. 

The two continuous-wave (CW) coherent Doppler lidar sys- 
tems (operating at wavelengths of 9.11 pm and 10.59 /xm) 
were developed by NASA/Marshall Space Flight Center 
(MSFC) [Rothermel et al ., 1995b]). Both beams were focused 
at about 50 m from the aircraft right-hand side with a typical 
sample volume of 600 m 3 (Table 1). Detection sensistivity for 
the 9.1 1-jxm and 10.59-^m lidars was ~6 x 10' 12 m' 1 sr' 1 and 2 
x 10' n m* 1 sr 1 , respectively. These sensistivity levels are 
mean values derived from an ensemble of signal-to-noise ratio 
data encompassing all GLOBE II flights. Thus the actual detec- 
tion level within a particular flight may vary around these val- 
ues. The lidar output signals were processed using two different 
data systems as described by Rothermel et al ., [1995b]. Both 
CW lidars also have been rigorously calibrated [Jarzembski et 
al ., 1995] using two different techniques that give internal 
consistency. 

The pulsed C0 2 coherent lidar system developed by NASA 
Jet Propulsion Laboratory (JPL) [Menzies and Trait , 1994] 
used a single wavelength of 9.25 jxm. Beam direction was 
selectable between either zenith or nadir positions. Pulse rate 
was 4 Hz. Data were stored on a pulse-by-pulse basis at a rate 


corresponding to 15-m range resolution, although the pulse 
length degraded the effective lidar range resolution to about 
100 m. These data were later edited, integrated over 100 s, and 
smoothed vertically over 150-m range increments. Detection 
sensitivity was «10' ! ! m' 1 sr' 1 at 1 to 2 km from the aircraft. 

The second pulsed lidar system (developed by NASA 
Goddard Space Flight Center (GSFC) [Spinhirne et al ., 1991] 
was based on a NdrYAG laser operating at a fundamental wave- 
length of 1.06 /xm. Raman shifting of the fundamental wave- 
length provided the 1.54-/xm transmitted pulse, while fre- 
quency doubling of the same wavelength gave the 0.532-/xm 
pulse. Beam direction was selectable between either zenith or 
nadir positions. Pulse rate used was 50 Hz and a range resolu- 
tion of 75 m selected. Aerosol backscatter threshold sensistiv- 
ity at 1.06 /xm and 0.53 mm was 1 x 10' 9 m' 1 sr' 1 and 5 x 
10' 9 m' 1 sr' 1 , respectively. Since Rayleigh backscatter can be 
of the order of or greater than the aerosol backscatter, this 
component was removed. Temperature profiles were used from 
the gridded National Meteorological Center model output 
closest to the lidar observation time and location [Spinhirne 
et al ., 1991]. These profiles were considered the best available 
at the time. 

3. Data Reduction and Analysis Procedures 

This section deals primarily with the preprocessing of the 
FSSP size distribution data, e.g., examining data quality, 
recalibrating the size bin edges, accommodating poor count 
statistics, and smoothing. To handle large quantities of FSSP 
and LOPC size distribution data efficiently, data visualization 
techniques were applied on a flight-by-flight basis. The com- 
puter software used was the JPL LinkWinds data visualization 
and analysis program [Jacobson et al ., 1994], which can dis- 
play data, after being suitably formatted, in a time-size cross- 
section plot with a false color display of the parameter of 
interest. Sample images of aerosol number in each size inter- 
val as a function of time are shown in Plate la (FSSP) and Plate 
lb (LOPC) for the Japan to Hawaii transit flight 14. Data 
displays in this format clearly showed any anomalous instru- 
ment responses. Any such problems were corrected and reveri- 



Plate 2. Contour plot of frequency of occurrence (percent) of FSSP-derived 9.1 l-/xm aerosol backscatter as a 
function of FSSP dn!d log D p averaging time for flight 14. 
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fied by comparing the two sets of data. These displays were 
also used to identify regions with unusual aerosol properties 
by verifing that both instruments observed similar features. 
One exception to this was when the aircraft flew through 
cloud. The FSSP continued to make measurements, while the 
LOPC was usually turned off. This can be seen in Plate 1 prior 
to 0000 and around 0200 UTC. 

3.1. LOPC 

Raw LOPC data were converted to concentration as a func- 
tion ot size for each channel. Corrections were applied for 
losses of large particles in the inlet plumbing and for gains of 
large particles in the virtual impactor (e.g., net gains of about 
a factor of 10 at 3 ^m diameter [Porter et al ., 1992]). The cor- 
rected distributions were smoothed over several size bins to 
reduce counting noise. Differences in the size distributions 
among the various preconditioning channels were used to 
derive size -distributed concentrations for sulfuric acid, ammo- 
nium sulfate/bisulfate, and refractory materials assuming that 
these constituents were externally mixed [ Clarke , 1991, 1993; 
V. Srivastava et al., manuscript in preparation, 1995). Data 
from the VI channel were used wherever possible in these 
computations because of their improved count statistics for 
larger particles. Results integrated over all particle sizes typi- 
cally showed a molar ratio of ammonium to sulfate near 0.3 
[Clarke, 1993). This value indicates a sulfuric acid/ammonium 
bisulfate mixture, referred to hereinafter as acid sulfate. 

The measured size distributions were then corrected for con- 
densed water that had evaporated between the ambient and the 
LOPC optical sensor, since the relative humidity in the sample 
volume was less than 10%. Ambient relative humidities in the 
troposphere (derived at flight level from a frost point hygrom- 
eter) generally varied from 10 to 70%. Water amounts were 
estimated using growth rules for acid sulfate particles [Tang, 
1980], assuming growth under nonfreezing conditions. 
(Observational evidence [Pueschel et al. , 1994] suggested that 
sulfate particles are probably supersaturated droplets at the 
temperatures encountered in these flights and hence remain 
unfrozen). Particle growth amounted to a factor of « 1.1 at 70% 
relative humidity. Estimated uncertainty in D/D 0 over mea- 
sured relative humidity range is of the order of ±0.05. 
Refractive index data for acid sulfate (which was set by the 
constant molar ratio) and water were combined using standard 
volumetric mixing rules. This allowed the real and imaginary 
aerosol refractive index to be derived at 0.55-, L06-, 9.1 1-, 
and 10.59-^m wavelengths. Since 9.11- and 9.25-/xm wave- 
lengths are very close and refractive index variations between 
them generally small, no 9.25-/im values were calculated. The 
derived refractive index values for each LOPC sample period 
were used to recalibrate the FSSP size distributions for that 
same period and to support the recalibrated FSSP data in 
backscatter computations at 1.06, 9.11, and 10.59 gm 
(section 3.2). Estimates of the aerosol growth and refractive 
index were not used in regions where the FSSP or the CW lidars 
indicated cloud, which would interfere with the LOPC data 
acquisition, or in highly variable aerosol conditions, which 
would interfere with the LOPC data processing. 

Doing a propagation of error calculation for ammonium sul- 
fate composition gives a square root of the sum (rss) of all 
uncertainties in derived visible backscatter of 15% at RH « 
15% and 25% at RH = 70%. In some instances, sampling 
statistics could increase the rss value. Furthermore, the use of a 
constant molar ratio and a constant growth factor across the 


LOPC size range could also lead to an error increase. For 
example, growth can be overestimated for large refractory par- 
ticles (dust) which are treated as if they were completely solu- 
ble. These errors will likely be greatest in high relative humid- 
ity because of the substantial aerosol growth under these 
conditions. Further work is being done to improve the 
physicochemical model by determining the molar ratio and 
refractive index independently for each size bin (V. Srivastava 
et al., manuscript in preparation, 1995). The size distributions 
and refractive index inputs from the improved LOPC aerosol 
model were used in this paper to calculate aerosol backscatter 
at the same lidar wavelengths used in the FSSP modeling 
study. 

3.2. FSSP 

The original FSSP size distribution data were measured with 
the probe calibrated using latex spheres (refractive index of 
1.585-0i at the 0.63-//m He:Ne laser wavelength). This value 
of refractive index is higher than that usually found for 
aerosols composed of sulfuric acid, ammonium sulfate, or dust 
whose real parts lie in the range 1.4 to 1.5. The original FSSP 
sensor calibration would then lead to size distributions with 
generally lower sizes and artifacts around the resonant region 
(0.8- 1.0 gm diameter). Simulations of a similar sensor, the 
PMS active scattering aerosol spectrometer probe (ASASP- 
100X), were shown to greatly distort the returned size distribu- 
tion if the correct refractive index is not used [Kim, 1995]. The 
distortion was greatest in the resonant region of the instru- 
ment’s response corresponding to 0.4- to 0.8-/im particle 
diameter. Strapp et al ., [1992] also reported that calibration 
changes need to be taken into account for the FSSP-300 when 
hygroscopic aerosols are measured in high relative humidities. 
Therefore the FSSP size bin limits were recalibrated using 
refractive index data (for 0.55 /im) derived from thermally seg- 
regated LOPC size distribution data (see section 3.1). The 
FSSP calibration procedure is based on Pinnick et al. [1981] 
where Lorenz-Mie calculations of the instrument’s response 
were done using the new refractive index data for each flight. 
(Note that the differences in refractive index for wavelengths 
0.55 and 0.63 gm are considered very small.) Typically, a 
single refractive index value is derived for each LOPC sam- 
pling period in a Bight, and so the FSSP size bin edges are cor- 
rected using a constant refractive index over each of these 
periods. Bin edge corrections amounted to less than 0.1 pm in 
the submicron size range up to 0.5 pm above this size range. 

Once the 10 s basic FSSP size distribution data are recali- 
brated, “boxcar” averaging is performed. The choice of averag- 
ing period needs to be carefully considered. Patterson and 
Bowdle [1991] examined the effect of averaging aerosol 
backscatter coefficient derived from the small sampling vol- 
ume of an earlier FSSP unit. They found from log-probability 
plots that as the averaging period went from their basic Lmin 
period to 10-min time resolution, the importance of the small 
tail in the backscatter distribution diminished and the 
backscatter became more lognormally distributed. That part of 
the cumulative probability curve above 50% showed little 
change. The major changes in the lower part of the curve 
occurred when going from l-min to 5-min averages. In this 
analysis a slightly different approach is taken. For each 
GLOBE II Bight a contoured frequency of occurrence (percent) 
plot of aerosol backscatter coefficient (calculated at 0.1 log 
P intervals) versus averaging interval at 9.11-j/m wavelength 
is made. An example is shown in Plate 2 for the Japan to 
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Hawaii transit flight (Bight 14). This contour plot indicates 
that for averaging periods between 20 s and 50 s the frequency 
of occurrence of /I values shows a strong tail below 1.5 x 
10' 12 m _l sr' 1 which recedes after the averaging period exceeds 
50 s. The presence of this small tail initially indicates an arti- 
fact arising from the use of too small of a sampling volume 
where insufficient particles in certain size bins are not sam- 
pled adequately. This is consistent with the results presented 
by Patterson and Bowdle after applying a wavelength conver- 
sion of 0.33 to obtain 10.6-/xm backscatter data. Averaging 
periods above 100 s reveal that the shape of the mode at high 
backscatter is mostly independent of averaging period. 
Similar trends in contour plots from several other GLOBE II 
flights suggest an optimal averaging period of 100 s. 

The next step in the data processing is to fit a monomodal 
lognormal curve to the 100 s averaged dn/d log (D p ) data. In 
the troposphere, particles formed from smaller particles by 
coagulation or by condensation of vapors, involve statistical 
fluctuations, and the resulting size distribution can be 
described by a lognormal distribution [Horvath et aL, 1990], 
Approximating the measured data with a modeled curve has two 
advantages. First, large statistical fluctuations will be 
appropriately weighted in bins with low counts (0, 1, and 2, 
...) at the larger size ranges (> 1 /xm). This approach reduces 
some of the large high-frequency fluctuations that would 
otherwise characterize backscatter estimates from the FSSP 
data. Second, the effect of multivalued response (artificial 
bumps and knees in the data [Pinnick et aL , 1981] occurring in 
the 1- to 2-/xm and 4- to 8-jxm diameter regions) can be 
minimized without resorting to bin amalgamation and 
subsequent loss of size resolution. The curve fitting is 
implemented with a least squares technique whose quality can 
be assessed by comparing the goodness of fit for the 
normalized number distribution and the normalized volume 
distribution and intercomparing the two as reported by 
Horvath et aL [1990]. The quality of the final fitting is also 
indicated by the correlation coefficient, which in many cases 
is above 0.98. 

A number of FSSP dvld log (D p ) plots from several Bights 
were examined and found to exhibit a single-mode structure. 
Therefore it was assumed in this analysis that the size distribu- 
tions above 0.35 /xm diameter were single mode since most of 
the measurements were made in the cloud-free troposphere. For 
the dust regime sampled near Tokyo, Japan, the monomodal 
shape of the size distributions was not significantly degraded 
where a single lognormal curve would not be appropriate. 
However, when cloud was encountered at Bight level, the curve 
fitting did not perform very well; the size data were Bagged so 
that no backscatter calculations were done with these data. 
(Also backscatter calculations for clouds required a different 
refractive index which the LOPC-derived indices may not 
reflect, because it was generally switched off when flying 
through cloud.) Geometric mean diameter for number distribu- 
tions generally fell between 0.4 and 0.7 /xm and the geometric 
standard deviation between 1.4 and 1.8. Size distributions 
with geometric mean diameter < 0.1 /xm and > 1.5 /xm and cor- 
relation coefficients equal to zero or negative were also 
Bagged and omitted from the backscatter calculation. It should 
be noted that the mode diameters derived from the FSSP data 
were significantly larger than those derived from the LOPC 
data, due to differences in the small-particle detection thresh- 
olds for the two instruments. Pueschel et aL [1994] summa- 
rizes the effects of small and large ammonium sulfate particle 


contributions to the aerosol backscatter coefficient, as a func- 
tion of wavelength. (Using ammonium sulfate instead of an 
acid sulfate composition would not significantly change the 
small-particle contribution.) For infrared wavelengths, con- 
tribution from the small diameter size range (< 0.35 /xm) 
amounts to less than 1%. However, for the 1.06- and 0.53-/xm 
wavelengths the aerosol backscatter coefficient can be under- 
estimated by as much as 22% and 33%, respectively. 

Finally, the calculated 100-s aerosol backscatter coefficient 
data are subjected to a running average using ±1 data points 
either side of the 100-s time point. It is found (see section 4) 
that there are many times when bin counts were not recorded 
above 1 /xm diameter during periods of low aerosol numbers 
and, consequently, the aerosol backscatter coefficient Buctu- 
ated widely even after the size data were fitted with an analyti- 
cal curve. Using a “boxcar” average with a 300-s period did not 
smooth out the large backscatter Buctuations where the corre- 
sponding coregistered CW lidar backscatter data indicated min- 
imal Buctuation. By doing a running average on the FSSP- 
derived aerosol backscatter. the larger fluctuations were 
considerably reduced with little loss of physical features. 


4. Results and Discussion 

4.1. Aerosol Number Comparison 

Statistical theory requires that FSSP particle counting obey 
Poisson statistics if the particle count from a given size bin is 
to be representative of the aerosol sample. FSSP count data 
from several Bights were examined during periods of low 
infrared backscatter and were tested against a null hypothesis 
that counts in a particular size bin obey Poisson statistics. In 
most cases the counts were found to be not significant at the 
5% confidence level for size bins corresponding to 0.8 /xm or 
above. Below this size the particle counts were generally sig- 
nificant giving the FSSP good count statistics. Therefore to 
independently assess this finding, particle number concentra- 
tions from the FSSP and LOPC size distribution data were 
compared for particle diameters greater than 0.8 /xm. To make 
this comparision, the FSSP dn/d log ( D p ) data are averaged 
over a LOPC sampling period. During a LOPC in-flight sam- 
pling period, each of the five channels was selected for about a 
minute resulting in a particular channel having a total sam- 
pling time of 2-3 min over a 10- to 15-min LOPC period. 
Hence, the averaging period for the FSSP data is about 5 times 
longer than an individual LOPC channel period. This approach 
was taken because the start and end times for the individual 
channels were not available during the LOPC sampling period. 
The LOPC volume sampling rate was about 1.7 times larger 
than the effective FSSP sampling rate and, with the virtual 
impactor (channel 4) in line, provided an additional gain of 
about 3 and 9 at size bins l and 3 /xm, respectively, for 12 km 
altitude [Porter et aL , 1992]. Considering all of these effects 
(including plumbing losses), the LOPC VI channel provides an 
effective sample volume of «16 times larger than the FSSP for 
the larger particles (=3 /xm). Note that this factor will increase 
at lower altitudes. The disparity in sample volumes was even 
greater when backscatter values were computed from 100-s 
FSSP samples instead of the 10- to 15-min LOPC samples. 

It was essential that these comparisons be limited to near- 
homogeneous aerosol regimes, so the LOPC and FSSP data 
were selected for comparison where CW lidar backscatter vari- 
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ance was minimal and no cloud encountered. Figure 2 shows 
examples of LOPC (VI) versus FSSP dn plots from four suit- 
ably homogeneous regions on three GLOBE II transit flights 
(flights 6, 8, and 14): the dn for each sensor are derived from 
integrating the respective averaged dn/d log { D p ) data over the 
diameter ranges 0.8- 1.0 /mi, 1-2 /xm, and 2-3 /xm. The maxi- 
mum number of plotted data points for a given size range is 
determined by the number of LOPC sampling periods occurring 
in the three selected time periods, e.g., four in flight 6 and five 
in flight 14. Generic estimates were made of FSSP dn values 
that corresponded to /V = 1, 10, and 100 total particle counts 
using the relationship dnld log ( D p )*A log (D p ) = NIK At. ( K is 
the FSSP effective volume sampling rate which varies with 
aircraft speed, and At is the sampling period.) Figure 2 shows 
three scales (one for each size group) where each decade count 
level occurs with respect to the dn scale. Binomial count 
statistics imply that for a 50% probability of getting a count 
and a required minimum accuracy of ±33%, the total counts in 
the sampling period need to be above 10. Also the case of 
small number of counts, the binomial distribution would sug- 
gest the likelihood of undercounting occurring rather than 
overcounting. 

Data for the smallest size grouping (0. 8-1.0 /xm) in Figure 2 
indicate a nominal 1:1 comparison. These data are mostly 
associated with high FSSP counts. Even below 10 counts there 
is only a weak bias toward the FSSP undercounting. For the I- 
to 2-/xm diameter range, dn data from all three flights lie above 
the 1:1 line showing a clear bias to FSSP undercounting. There 
is a suggestion that the bias increases as the count rate 
decreases toward unity, with bias factor ranging from five to 
eight. The 2- to 3-/xm diameter range showed mixed results; 
five of the compared dn data show a bias. In this size range, 
relative humidity corrections, or correction for the large parti- 
cle loss to the LOPC particle size, may introduce additional 
errors masking the FSSP undercounting. 

The fact that the FSSP sensor is underestimating dn/d log 
(D p ) for particles above 1 /xm diameter is discussed further in 
the following section on infrared aerosol backscatter data 
comparisons. It is this size range that contributes greatest to 



Figure 2. Cross plot of LOPC (with virtual impactor) dn 
versus FSSP dn data for three selected size ranges taken from 
flight 6, 0036-0154 UTC, flight 8, 0129-0208 UTC, and 
flight 14, 0300-0345 UTC and 0415-0445 UTC. Values of dn 
corresponding to FSSP decade total particle counts are 
indicated at the top of the Figure by the FSSP count scale for 
each size range. 


the differential aerosol backscatter coefficient at these wave- 
lengths [Srivastava et al., 1992]. 

4.2. Aerosol Backscatter Coefficient 
Comparisons: Case Studies 

Results presented here are based on transit flights 6 and 14. 
Flight 6 was chosen as a representative of a meridional flight 
across the equatorial region, while flight 14 represents a zonal 
flight in the northern hemisphere. 

Occasional breaks in the FSSP-derived shown in the time 
series plots for the two transit flights are the result of criteria 
applied to the lognormal fitting of the size distribution data, 
as discussed earlier. For example, no fitting is done when 
cloud was encountered because the appearance of a cloud parti- 
cle mode made fitting a single-mode lognormal curve inappro- 
priate: this study was confined to clear air aerosol backscatter. 
Since the 9.11- and 9.25-/xm wavelengths are very close and 
refractive index variations between them were usually small, 
no backscatter computations were performed at 9.25 /xm; the 
9. 1 l-/xm FSSP-derived backscatter data, /^fssp. 9. id* were used 
to compare with the 9.25 /xm pulsed JPL lidar backscatter, 
Ajpl, 9.25). 1° case °f ^e LOPC-derived fi calculations for 
0.53-, 1.06-, 9.1 1-, and 10.59-/xm wavelengths the refractive 
index values were calculated using a physicochemical mixture 
model which determines the molar ratio for each size interval 
rather than a constant molar ratio across the whole size 
distribution. Hence, a growth factor was inferred that is more 
appropriate for that size and composition (V. Srivastava et al., 
manuscript in preparation, 1995). Aerosol backscatter data for 
several wavelengths from each of the two transit flights are 
now discussed. 

4.3. Hawaii to American Samoa Transit Flight 

Flight details. Flight 6 covered the leg from Honolulu, 
Hawaii, on a track south-southwest to Pago Pago, American 
Samoa (see Figure 1), on May 20-21, 1990, from 2105 to 
0208 UTC with a constant pressure altitude maintained at 
approximately 8.5 km for the duration of the flight. At about 
2130 UTC the aircraft encountered westerly winds which 
reached a maximum of 20 m s* 1 at 2230 UTC, while at 2245 
UTC, enhanced cloudiness and some convection was noted at 
flight level. This cloudiness was associated with a band of 
cloud to the east and northeast which could have had frontal 
characteristics. At 2300 UTC a 400-km band of light winds 
and deep convection associated with the Intertropical 
Convergence Zone (ITCZ) was encountered. South of the ITCZ 
at 0020 UTC some isolated convection was observed as evi- 
denced by the high CW lidar backscatter and high relative 
humidity values at aircraft level. After this period (0100 UTC), 
winds become southerly at less than 5 m s' 1 which generally 
prevailed until the end of the flight. Equator crossing occurred 
at 0000 UTC. Figures 3a-3c show the aerosol backscatter 
coefficient computed from the FSSP size data for the wave- 
lengths 1.06, 9.11, and 9.25 /xm, respectively. Very little 
GSFC lidar data were available for a meaningful comparison at 
0.53 /xm, and the 10.59-/xm CW lidar data were very noisy; 
consequently these data were not used. 

1.06 /xm. 0(FSSP, 1.06) an d 0(GSFC. 1 . 06 ) data are shown in 
Figure 3a. Over the time period 2120-2300 UTC, temporal and 
absolute backscatter trends are followed reasonably well even 
below the GSFC lidar threshold line. The two peaks seen in the 
Agsfc. 1 . 06 ) data at 2145 and 2200 UTC were the result of the 
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Figure 3. Temporal variation of measured and modeled 
aerosol backscatter coefficient (£) during flight 6 for 
wavelengths (a) 1.06 gm, (b) 9.11 gm, and (c) 9.25 /im . 
Dashed line in each panel represents the lidar detection 
threshold. 

top portions of an aerosol plume (seen in the time-altitude 
backscatter cross-section image) entering the averaging 
region below the aircraft. Backscatter values within these 
plumes were about 2-3 x 10* 8 m' 1 sr* 1 . This plume originated 
from the marine boundary layer (MBL) near the ITCZ (latitude 
12°N) and extended up to about 7 km. ^(gsfc, i. 06) data around 
2220 to 2240 were heavily edited to remove cloud returns and 
are below the threshold line which resulted in some scatter in 
the data. 


For the period 2340-0030 UTC the /3 ( gsfc, i. 06) data are 
scattered to either side of the FSSP-derived data. The time- 
altitude j3(GSFC, 1.06) cross-section image indicates a region of 
aerosol inhomogeneity and editing of cloud returns. In the 
third comparison period (0110-0145 UTC), j8 (G sfc. i.06) data 
show a similar temporal trend to the £ ( fssp. i .06) data b ut are 
about a factor of 5 higher. In this case the nadir GSFC data 
were available for averaging over only the first 0.25 km from 
1 km below the aircraft instead of the usual 1-km depth. 
Having a smaller number of points in the averaging may have 
contributed to a small part of this difference, although the 
points included were nearest to the flight level. The time- 
altitude ^gsfc, 1 . 06 ) cross-section image for this period does 
indicate a very uniform aerosol backscatter. However, this 
does not rule out the existence of spatial aerosol inhomogene- 
ity between flight level and 1 km below. The LOPC differential 
backscatter plot (not shown) for particle diameter range 0.16- 
0.35 indicated only 25% contribution to the total 1.06-/mi 
backscatter. On the other hand, the condensation nuclei (CN) 
plot (also not shown) also indicated an increase in the ambient 
CN from 400 to 1200 cm' 3 with very little refractory CN pre- 
sent, but at the most, this would imply only an additional 
small increase in differential assuming no aerosol composi- 
tion changes. The ^ (L opc, i.06) data showed best agreement 
here when an acid sulfate composition is selected. 

9.11 gm. Figure 3b shows the aerosol backscatter coeffi- 
cients derived from FSSP size distribution data for 9. 11-gtm 
wavelength using refractive index data derived from the LOPC 
thermally segregated size distribution data, along with the cor- 
responding CW lidar and LOPC-derived aerosol backscatter 
data. In the early period from 2115 UTC until cirrus cloud was 
encountered at 2250 UTC, the $ F SSP. 9.1 1) data closely follow 
the /I(cw, 9.H) data except in the period 2200 to 2220 UTC. 
The underestimate in £ ( fssp. 9.1 d data for these low values 
occurred consistently in these analyses as well as in other 
analyses not reported here. After the cirrus cloud region the 
/ 2 (cw, 9 . 11 ) data lay above the /^fssp. 9.1 1 > data *° r the rest °f 
the Bight except at 0050 and 0120 and showed much less tem- 
poral variation. The jS (LO pc, 9 . 1 1 > data followed much closer to 
the /^cw. 9 . 11 ) data m absolute levels and also exhibited less 
temporal variation than the $fssp, 9. 1 1 > data. One of the rea- 
sons for the differences seen here is the count statistics for 
FSSP size data above 1 //m diameter. Particles in the 0.8- to 3- 
/im-size range dominate the backscatter in these infrared wave- 
lengths [Srivastava el aL, 1992], and so there will be the ten- 
dency to underestimate the computed backscatter if these parti- 
cles are not adequately sampled. Using FSSP averaging times 
greater than 100 s would not improve this situation: the FSSP 
would have to actually sample for much longer than 100 s to 
improve the large particle count statistics, or the sampling 
procedure would have to be modified to concentrate these par- 
ticles into the FSSP sampling volume. On the other hand, the 
LOPC, already equipped with virtual impactor and using sam- 
pling times of several minutes, had better counting statistics 
(see section 4.1). Agreement between ^ cw . 9.1 d and 
A LOPC, 9.1D was ver y g°°d m l° w backscatter conditions over 
the period 2335 to 0200 UTC. Except for the cirrus cloud 
period the aerosol composition derived from the LOPC data 
analysis is best represented by the acid sulfate composition 
over all size ranges (V. Srivastava et aL, manuscript in 
preparation, 1995). The excursions in FSSP-derived at 0050 
and 0120 UTC can be explained as follows: If the composition 
of large particles (>1 fim) over this period does not change, 
then the lidar and LOPC backscatter is expected to not show 
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much variation. On the other hand, if the submicron particle 
composition changes, then the FSSP-derived backscatter could 
vary significantly because FSSP is only sampling this range 
of particles adequately and not the larger particles. For 
example, the submicron particles at 0050 and 0120 UTC may 
possess more ammonium sulfate which would result in a much 
higher backscatter, because of the higher imaginary 
component of the refractive index. Analysis of other GLOBE II 
flights revealed similar results (D. R. Cutten et. al., 
manuscript in preparation, 1995). 

A value for the underestimated differential aerosol backscat- 
ter (df$/d log (D p ))* A log D p over a given size range has been 
derived. On examining the dn data for the 1 - to 2-/xm size range 
(Figure 2), a deficiency factor of about 6 was estimated for the 
FSSP data. An average value of dpid log ( D p )*A log D p for a 
weak acid sulfate aerosol is about 5 x 10 " 12 nr 1 sr * 1 and when 
multiplied by this deficiency factor, becomes 3 x 10 " 11 m " 1 
sr' 1 . From Figure 3b the average J3 (fssp, 9 . it) value for the time 
period 0030-0155 UTC (excluding the two major peaks) was 
about 1 x 10 ' 1 1 m " 1 sr* 1 . Adding in the dp/d log (D p )*A log D p 
value gives a new ^(fssp, 9 .ii) value of 4 x 10- 1 ! m~ 1 sr* 1 , 
which is much closer to the £<cw, 9 . 1 1 ) data. 

9.25 jun. Figure 3c shows the comparison of /3 ( j PL 9.25) 
data and ^(fssp. 9.1 1) data which indicated varied agreement 
along the flight. In the earlier part of the flight (2130-2250 
UTC) the 9 25) data were generally above the / 3 (fssp, 9.1 1) 
data with an ill-defined temporal trend. Cirrus cloud which 
extended down to flight level was responsible for the high 
backscatter in the period 2250 to 2340 UTC. For the time 
period 2350 to 0200 UTC again the /?(fssp, 9.1 d data were gen- 
erally below the £(jpl, 9.25) data. The time-altitude ^ JPL 9.25) 
cross-section image indicated fairly uniform backscatter 
within the averaging region with several very low backscatter 
areas around 0030 UTC. Apart from the low P regions at 0030 
the p( j PLi 9 25) data follow the /?<cw, 9.1 d data. For the same 
reason discussed above, poor count statistics for large 
particles are mostly responsible for the FSSP underestimating 
the backscatter, whereas the 0(lopc, 9.10 data, however, are in 
closer agreement with the $jpl, 9.25) data. 

4.4. Japan to Hawaii Transit Flight 

Flight details. Flight 14 covered the leg from Tokyo, 
Japan, to Honolulu, Hawaii (see Figure 1 ), on June 3-4, 1990, 
from 2255 to 0530 UTC with a constant pressure altitude 
maintained at 8.5 km for the duration of the flight. Detailed 
meteorological conditions encountered along the flight are 
given by E. W. McCaul et al. (manuscript in preparation, 
1995) and only a brief synopsis is given here. The flight 
began near a midlatitude jet located south of Japan, and after 
several hours into the flight, several shortwave troughs were 
traversed. During this period, cirrus cloud was generally 
located above the aircraft but was at flight level during the 
periods 2300-0000 UTC and 0150-0215 UTC. On the 
approach to Hawaii, subtropical high conditions prevailed 
with little cloud present at or near flight level. Figures 4a-4e 
show the aerosol backscatter coefficient computed from the 
FSSP size data for the wavelengths 0.53, 1.06, 9.11, 9.25, 
and 10.59 /xm, respectively, along with the lidar and LOPC- 
derived /Jdata. No lognormal fitting was done to the FSSP size 
data for the periods cirrus cloud was traversed. Several other 
small gaps are seen in the plot where the fitting could not be 
done reliably. 

First, the discussion will briefly dwell on the composition 
of aerosol particles encountered during this flight. A more 



Time (UTC) 



Figure 4. Temporal variation of measured and modeled 
aerosol backscatter coefficient ( P ) during flight 14 for 
wavelengths (a) 0.53 /xm, (b) 1.06 /xm, (c) 9.1 1 /xm, (d) 9.25 
/xm, and (e) 10.59 /xm. Dashed line in each panel represents the 
lidar detection threshold. 


detailed discussion on the modeling of the aerosol composi- 
tion encountered during this flight, using infrared aerosol 
backscatter data to constrain the physicochemistry model, 
will be published (V. Srivastava et al., manuscript in prepara- 
tion, 1995). Both measured and derived aerosol backscatter for 
9.1 l-/xm and 10.59-/xm wavelengths (see Figures 4c and 4e) 
indicated that for the initial portion of the flight up to about 
0230 UTC (excluding the period of cirrus cloud), a partially 
ammoniated continental dust aerosol represented the particle 
regime. After 0230 UTC, varying amounts of ammonium sul- 
fate, characteristic of an advected continental plume, were pre- 
sent, diluted with clean tropospheric air. During 0350-0515 
UTC the aerosol was dominated by sulfuric acid with some par- 
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Figure 4. (continued) 


tial ammoniation comfirmed by the 9.11-/zm and 10.59-/zm 
CW lidar p data for this period as having a ratio near 2. During 
the last period the composition is consistent with an aerosol 
type generally associated with low backscatter in the free tro- 
posphere of a subtropical high [Srivastava et ai, 1995], 

0.53 and 1.06 gm. The comparison of FSSP-derived p 
data with GSFC P data at these two wavelengths showed rea- 
sonable agreement over two portions of the flight, while sev- 
eral other regions showed large disagreement. One of the 
largest discrepancies occurred during 0050-0120 UTC where p 
estimates from both the FSSP and the LOPC were initially 
much smaller than the GSFC data, especially at 0.53 /zm, but 
gradually approached the GSFC data at both 0.53 and 1.06 /zm 
later in the period. Flight level volatility data initially showed 
a predominately refractory sootlike CN aerosol with little 
change in CN concentrations until later in the period. During 
this same period, LOPC data showed a systematic increase in 
larger particles most likely composed of mixed acid sulfates. 
However, up to 30% of the LOPC-derived P values was con- 
tributed by particles between the LOPC threshold (0.15 /zm) 
and the FSSP threshold (0.35 /zm). GSFC data initially showed 
much higher /3at 0.53 /zm than at 1.06 /zm. However, measured 
P at 1.06 /z m increased steadily, while p at 0.53 /zm decreased 
steadily. All of these trends are consistent with a systematic 
change in aerosol type from smaller particles early in the 
period to larger particles later in the period, both at and below 
night level. Discrepancies between FSSP and GSFC p values 
during this period were probably due to high concentrations of 


small particles below the FSSP detection threshold, combined 
with vertical gradients in a stratified plume such as might be 
expected in aged urban pollution downwind of Tokyo. 

In the period 0145-0220 UTC the aircraft passed through 
some cloud as was evident from the FSSP size distribution 
image and CW lidar data plots. However, below the aircraft 
(>1 km) there was no evidence of returns from cloud in the 
pulsed lidar backscatter data. Therefore comparisons of the 
FSSP data with the GSFC lidar data are not meaningful because 
they were sampling different scattering media. The LOPC- 
derived p show some agreement with the pulsed lidar data 
probably because the cloud droplets impacted on the walls of 
the plumbing or evaporated before they reached the LOPC 
sample volume. 

For the last period (0300 to 0515 UTC) the 1 .06-/zm com- 
parison indicated reasonable agreement between GSFC- and 
FSSP-derived p data except for the fall in p occurring at 0330 
UTC. The l.06-/zm GSFC time-altitude cross-section backscat- 
ter image data for this flight indicated a 300- to 500-m-thick 
aerosol layer present 1-2 km below the aircraft flight level and 
its upper interface inclined upward in the flight direction. As 
the aircraft progressed along the flight track, the nadir- 
viewing GSFC lidar was probing less and less of this layer 
resulting in a lower averaged p leveling out near the threshold 
line. Once the aircraft cleared this high-aerosol backscatter 
regime, the two p estimates merge in the low backscatter 
regime even though that the GSFC p data fall below the 
threshold level for this wavelength. The short comparison 
around 0500 UTC, still in this low-level backscatter regime, 
indicated that the GSFC and LOPC P data lie above the FSSP- 
derived p data before aircraft descent. This difference can be 
accounted for by the large contribution (>20%) to P from the 
0.16- to 0.35-/zm-diameter particles that the LOPC sensor was 
able to measure. The 0.53-/zm FSSP-derived and GSFC j3data 
comparison covered a much shorter period (0300 to 0330 
UTC) and showed agreement only in the latter part of the 
period before the large drop when entering the subtropical 
high region. LOPC differential backscatter data for this period 
indicated only a small contribution (< 1 5%) to p from the 
particles in the size range 0.16 to 0.35 /zm diameter. In both 
the 0200 and the 0300 hour periods, ratios of GSFC 0.53- and 
1.06-/zm backscatter data indicated that the wavelength 
exponent fluctuated between 1 and 2.3 which indicated the 
existence of large and moderate size particles, respectively. 

9.11 and 10.59 /zm. In these wavelength comparisons 
(Figures 4c and 4e) both the CW lidar and the FSSP-derived p 
were obtained at flight level, thus minimizing the effects of 
aerosol inhomogeneities. In this flight the large temporal 
trends are closely matched although smaller changes are not 
followed as well. Some of the smaller FSSP-derived P data fluc- 
tuations arise from count-sampling statistics in a small sam- 
pling volume, whereas the CW lidar does not necessarily 
exhibit such fluctuations because of its much larger sampling 
volume (a factor of ~6xl 0 5 larger). On longer timescales, 
measured and modeled 9.1 l-/im p values generally agree within 
a factor of 2 except for a factor of 3 discrepancy during the 
period 0415-0445 UTC. A similar comment can be made for 
the 10.59-/zm wavelength comparison (Figure 4e). Note that 
the period 0000-0100 UTC showed better absolute agreement 
at 10.59 /zm than 9.11 /zm although a different temporal trend 
was evident in the 10.59-/zm CW lidar p data. 

The biggest difference in the 9.11- and 10.59-/zm p data 
occurred at 0415-0445 UTC during the traversing of the sub- 
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tropical high on the approach to Hawaii. This can be 
explained as follows: The aerosol backscatter in this period 
was dominated by a low concentration of particles around I fi m 
diameter which the FSSP was unable to sample adequately but 
were sampled more accurately by the CW lidar beam because of 
the five order increase in sampling volume. The LOPC-derived 
P data are in good agreement with the ^ (CW 9 , and 
P( CW, io.59) data. The false color size-time cross-section plot 
of dnld log ( D p ) for this flight (see Plate la) confirms that the 
FSSP concentrations above 1 ^m diameter were almost 
neglible during this time period, while the LOPC 
concentrations were as high as 0.01 cm* 3 . A calculation of 
9.1 1 -fm\ differential backscatter coefficient, dP/d log (D p )*A 
log D p (as for flight 6 infrared P data), confirmed this finding, 
yielding a corrected P value around 2.3 x 10“ 1 1 m' 1 sr' 1 for this 
period. A deficiency factor of 2 in dn was derived from the 
FSSP-LOPC dn plot. Aerosol composition for this region 
indicated relatively more acidic sulfate aerosol present than in 
the case examined in flight 6, resulting in a slightly higher 
dPld log (D p )* A log D p value for the 1- to 2-pm size range (V. 
Srivastava et al., manuscript in preparation, 1995). 

9.25 fxm. Comparison of P( pssp, 9 . 1 1 ) and P(spl, 9.25) data 
(Figure 4d) showed much more variability than any of the 
other wavelengths compared in this Bight. Some of this vari- 
ability arose from aerosol inhomogeneity near the flight path 
because the j5 (JPL 9 25) data are measured either 1 km below or 
above the aircraft flight level. Lidar backscatter near 0000 
UTC was due to a high backscatter layer in the 6- to 7-km alti- 
tude region. In the time period 0130-0145 UTC the high P 
values are due to cirrus occasionally extending down nearly to 
flight altitude. Thus the situation at flight level was not 
always representative of what the JPL lidar was seeing. JPL 
trends here also tend to correlate with GSFC trends, thus sup- 
porting aerosol inhomogeneity. The closest agreement is for 
the period 0400-0500 UTC where p showed an overall trend to 
lower values. 

To summarize, the infrared wavelength-measured and mod- 
eled comparisons of aerosol backscatter near or above 10' 9 
m' J sr' 1 showed reasonable agreement. Below this p level, 
agreement is not so good. In most cases, poor count statistics 
for particles greater than about 1.0 /im diameter are contribut- 
ing to a lower and more fluctuating FSSP-derived p in low- 
aerosol backscatter regions. 

4.5. Aerosol Backscatter Coefficient 
Comparisons: Other Transit Flights 

In this section, three other transit flights are examined: two 
(flights 8 and 10) cover the southern hemisphere, while the 
third (flight 12) is another meridional flight across the equato- 
rial region. Data from these flights have been combined to 
provide cross plots of FSSP-derived and lidar-derived P data 
rather than the time-series plots shown in the previous sec- 
tion. LOPC-derived p are not included in these plots because of 
the much longer averaging time used to sample the aerosol and 
variations in CW lidar p would be smoothed out. Data from 
flights 6 and 14 have also been included in the plots for com- 
parison purposes. Any cloud returns from the CW lidars were 
excluded and all remaining data smoothed with a 300-s running 
average. Flight details are as follows: 

Flight 8. This flight flew southwestward from Papeete, 
Tahiti, to Christchurch, New Zealand, on May 23-24 begin- 
ning at 2258 and ending at 0524 UTC. Flight altitude varied 
between 8 and 12 km. Scattered convection gradually increased 


in coverage early in the flight until at 0040 UTC the aircraft 
was forced to climb to 12 km to fly over the storms. At about 
0230 UTC the aircraft descended to 9 km and traversed a region 
suggestive of a midtropospheric front (latitude 33° S). After 
descending further to 8 km near 0300 UTC, the aircraft entered 
into a region of anticyclonic flow. After 0300 UTC the aircraft 
began to overfly increasing amounts of lower-tropospheric 
convection and was forced at 0430 UTC to ascend to 12 km 
near latitude 44° S. Westerly winds of the order of 60 m s‘ l 
were encountered in a jet core at this altitude shortly before 
landing. 

Flight 10. This flight went from Christchurch, New 
Zealand, to Melbourne, Australia, on May 28 starting at 0022 
and ending at 0805 UTC. Altitude varied from about 7.5 km to 
8.5 km with the lowest level being at the extreme latitude of 
56 S at 0330 UTC. Low-level cloudiness was widespread in 
this flight along with some shallow convection. Near latitude 
55 S a 70 m s‘ l westerly jet was encountered before landing. 

Flight 12. This flight commenced at Darwin, Australia 
(0058 UTC), and terminated at Tokyo, Japan, on May 31 
(0723 UTC). At around 0340 UTC a clean «200-km-wide N-S 
midtropospheric air mass was encountered just north of the 
ITCZ. However, by 0430 UTC the air mass had changed over to 
a predominantly continental air mass containing Asian dust in 
a low humidity environment by which time a strong westerly 
upper level jet was encountered. The plume persisted until the 
descent into Tokyo. This dust feature was strongly evident in 
the time-altitude backscatter cross-section images for the 
GSFC and JPL pulsed lidars (E. W. McCauI et al., manuscript in 
preparation, 1995) which also includes a more detailed 
description ot the meteorological conditions encountered. 

9.11 and 10.59 gixn. The 9.11- and 10.59-^m data (Figures 
5a and 5b) at P values above ~2-3 x 10' 10 m' 1 sr’ 1 are clustered 
either side of the 1:1 line, which indicates that the count 
statistics for the FSSP appear adequate for these backscatter 
values. However, there is still a wide range of scatter in this 
backscatter region which is probably the result of different 
factors. As noted earlier, the disparity between the sampling 
volumes of the FSSP and CW lidar is quite large and probably 
highlights aerosol inhomogeneities resulting from refractive 
index changes on a small scale. Furthermore, the FSSP sam- 
pling volume, although located on the same side of the aircraft 
as the CW lidar sampling volume, did not fall within the lidar 
sampling volume. 

For p < 10' 10 m' 1 sr' 1 data points at the two wavelengths for 
all flights lie mostly above the 1:1 line (Figures 5a and 5b) 
with the most deviation away from this line at the very low 
levels, confirming the earlier conclusion that poor FSSP count 
statistics imply lower FSSP-derived backscatter values. Most 
data from flights 8 and 10 show aerosol backscatter lying 
within this region of low backscatter. The shaded band in the 
two figures indicates the range of FSSP-derived P when only 
one count is registered in one of the FSSP size channels 
between 0.7- and 1.3-/xm particle diameter; all other channels 
are assigned zero counts. Note that these band limits depend 
on the refractive index of the aerosol being sampled and the 
FSSP calibration constant. Therefore the shaded band repre- 
sents an “effective backscatter detection” threshold for the 
FSSP in regions with occasional large particles, roughly anal- 
ogous to the more rigorous detection thresholds for each of the 
lidars. The shaded area is representative for flight 14 and indi- 
cates that the scatter of backscatter at these low values can be 
mostly attributable to only a few counts being available. One 
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Figure 5. Cross plots of lidar-derived versus FSSP-derivcd 
aerosol backscatter coefficient (/3) at (a) 9.1 1 /im and (b) 10.59 
jim for Bights 6, 8, 10, 12, and 14. Horizontal dashed line in 
each panel represents lidar detection threshold, while diagonal 
dashed line pair lies at ±3 factor from the 1:1 line. Shaded area 
represents range of FSSP-derived tor one bin count in size 
bins between 0.7 and 1.3 (see text). 


should also note for Eights 8 and 10 that there are a few 9.1 1- 
and \0.59-fim backscatter data falling below the CW lidar 
threshold line and the t:l line, whereas for Eight 12 data, this 
is not the case. Possible reasons are as follows: First, any data 
below the threshold line should be considered noisy and have a 
large error [Rothermel et aL , 1995b]: Eight 10 data at 10.59 
jxni, relatively, having the most points below the threshold 
line. Second, those data points below the 1:1 line and above 
the CW lidar threshold line are probably due to a combination 
of spatial variability (since sensor sampling volumes differ by 
a factor of 6 x 10 5 (Table 1)) and sampling statistics. 


4.6. Free Troposphere Background Aerosol Mode 

Finally, all the 9 . 1 1 ) an d ^(cw, 9 . 1 1) data at Eight 

level and /3 (jpl 9 25) data within 1 - 2 km of the aircraft (for all 
five Eights) are replotted in the form of three histograms 
(Figures 6a-6c). Each plot represents data that has been further 
subdivided into northern and southern hemisphere categories, 
and Table 3 summarizes the low and high fi modes associated 



-12 -11 -10 -9 -8 -7 -6 


Log 9.1 1 pm FSSP backscatter (m 1 sr ^ 



-12 -11 -10 -9 -8 -7 -6 

Log 9.1 1 pm CW lidar backscatter (m 1 sr ') 



Log 9.25 pm JPL lidar backscatter (m 1 sr ') 


Figure 6. Frequency histograms of mid-infrared aerosol 
backscatter coefficient for the northern and southern 
hemispheres derived from (a) FSSP, (b) continuous wave (CW) 
lidar, and (c) Jet Propulsion Laboratory pulsed lidar sensors, 
respectively. The FSSP and CW lidar data are for Eight level; 
JPL data are within 1-2 km above and below Eight level. 
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Table 3, List of Modes Seen in Free Troposphere Infrared Aerosol Backscatter Histograms 


Northern Northern Southern Southern 

Hemisphere Hemisphere Hemisphere Hemisphere 

Sensor (Low p Mode) (High p Mode) (Low p Mode) (High p Mode) 


GLOBE II 

FSSP-derived P 

1x10-' 1 

IxlO’ 9 

IxlO" 

- 

GLOBE II 

(9.1 1 /rm) 
CW lidar P 

3x10-" 

IxlO’ 9 

3x10" 

weak 

GLOBE II 

(9. 1 1 ji m) 
JPL lidar p 

3x10' 1 

8x10" 10 

4x10" 

weak 

GLOBE I 

(9.25 /im) 
FSSP-derived P 

IxlO' 1 1 


3x10" 



(~ 10 /im) 






with each histogram. Each histogram for the northern 
hemisphere clearly shows two modes with an intermodal gap 
near 1 x 10' 10 . The high-backscatter mode in all plots has a 
common center around 1 x 10' 9 m' 1 sr" 1 . However, the FSSP- 
derived low backscatter or background mode (1 x 10" 11 nr 1 
sr' 1 ) is clearly lower than the corresponding mode in the other 
two northern hemisphere histograms. A similar shift is also 
seen in the southern hemisphere histogram plots. The 
background modes (3 x 10' 11 nr 1 sr' 1 ) for aerosol backscatter 
measured with the CW Doppler and JPL lidars are in qualitative 
agreement with background modes observed in other CO 2 lidar 
aerosol backscatter data obtained in the free troposphere 
[Rothermel et al . , 1989]. Hence, the lower value for the back- 
ground mode observed for the FSSP-derived data must be the 
result of the FSSP sensor underestimating the 0.8- to 3-jum 
diameter particles as reported earlier. The factor of 3 is also the 
value found for regions of flights 6 and 14 where dpid log 
(D p )* A log (D p ) over this size range needed to be increased. A 
similar histogram of aerosol backscatter (~10-jim wavelength) 
has been presented by Pueschcl et al. [1994] derived from a 
number of randomly selected FSSP size distributions measured 
in GLOBE I. (These calculations were based on original Latex 
sphere calibration for the FSSP and an ammonium sulfate 
aerosol refractive index.) In that study the northern and south- 
ern hemisphere modes were located at 1 x 10' 11 m' 1 sr* 1 and 
3 x 10' 1 1 m' 1 sr' 1 , respectively: the latter mode is a factor of 3 
higher than reported here for the 1990 southern hemisphere 
data. Using the P correction factor of 3 estimated earlier, 
PuescheLs background modes for the northern and southern 
hemispheres become 3 xlO* 11 and 6 x 10* 11 m* 1 sr' 1 , respec- 
tively. However, since Pueschel ct al. ’s P calculations are not 
based on an acid sulfate aerosol and wavelength is higher than 
used here, their background mode would be enhanced by 
20-30%. The higher-backscatter mode seen in the southern 
hemisphere during GLOBE I is similar to that measured with 
the other lidars during this austral 1989 spring mission. The 
northern hemisphere FSSP-derived background modes, how- 
ever, are now the same. 

The upper mode at 1 x 10* 9 m' 1 sr’ 1 present in the northern 
hemisphere (Figures 6a-6c) is mostly the result of aerosol 
backscatter measured during flight 14 with a small contribu- 
tion from the dust region encountered in Bight 12. There was 
no contribution to this mode from flight 6. The corresponding 
£(LOPC, 9.11) histogram (not shown) covering the northern 
hemispheric region confirmed that sulfate aerosols, possibly 
mixed or coated with dust, contributed significantly to the 
backscatter in this mode. Note, however, that this mode was 
very weak in the southern hemisphere. 


5. Conclusions 

One of the significant conclusions from this work is that 
when the FSSP optical particle counter is used with short 
sampling periods (<10 s) in clean atmospheric conditions, it 
underestimated aerosol particle concentrations within the size 
range 0.8 to 3.0 jum diameter. Particle counts within a 
particular size bin should conform to Poisson statistics if they 
are to be representative of the aerosol, which was not the case 
for size bins above 0.8 gni diameter because of the small 
sample volume. Evidence for undercounting was confirmed by 
( 1 ) direct comparison of dn averaged over selected time periods 
and integrated over several size ranges for FSSP and LOPC 
(with virtual impactor) and (2) comparison of infrared aerosol 
backscatter estimate at Bight level measured with CW C0 2 
lidars and derived from the FSSP. Pfssp data are underestimated 
below 10' 10 m' 1 sr' 1 at wavelengths betw-een 9.11 and 10.59 
fim because this is the size range which contributes most to 
the infrared aerosol backscatter. Backscatter derived from the 
FSSP have been estimated to be as much as a factor of 3 too 
low when backscatter falls between 10* 11 nr 1 sr" 1 and 10' 10 m" 1 
sr' 1 . The large scatter seen in the Pfssp data can also 
attributed to insufficient counts in the 0.8- to 3.0-/im-diameter 
size range. On the other hand, LOPC (VI) provided better 
agreement in this backscatter range. This is because the LOPC 
had a net gain over the FSSP in this size range from the use of 
a virtual impactor and larger sampling volume. 

For p > 10' 10 nr 1 sr' 1 the agreement between the FSSP- 
derived and Bight level lidar data is better because FSSP count 
statistics have greater statistical significance, even though 
some of the data were noticeably scattered. The scatter likely 
arises from the large CW lidar sampling volume exhibiting 
aerosol microphysics variations (from small refractive index 
changes) within the volume that causes aerosol backscatter to 
Buctuatc. CW lidar aerosol backscatter close to I0" M nr 1 sr' 1 
showed better agreement with the FSSP-derived data because 
particles smaller than I jim contribute most to the infrared 
backscatter under these conditions and the FSSP had 
reasonable count statistics for these size channels. This 
conclusion about underestimating particles in the 0.8- to 3-/im 
size range means that to adequately obtain a representative 
size distribution one has to either increase the FSSP sampling 
time considerably or devise some means of concentrating 
these particles into the FSSP sampling volume to maintain 
short sampling times. One would still need to apply an 
analytically fitted lognormal curve to the data to reduce the 
inBuence of multivalued bin counts. 

The comparison of the FSSP-derived and GSFC lidar 0.53- 
and 1.06-jim aerosol backscatter indicated that the FSSP- 
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derived data had generally underestimated the lidar data as 
expected. The 0.35-^m size cutoff for the FSSP meant no 
particles were sampled adequately below this range, which 
contributes about 20% of the backscatter at 1 .06 pm and 33% 
at 0.53 pm [ Pusechel et al ., 1994]. Many differences between 
the GSFC lidar and the FSSP-derived backscatter were 
attributable to aerosol inhomogeneity as sensor sampling 
volumes were not colocated, as well as there being a large 
disparity between the two sample volumes. The P^fssp. 9.U) 
and P(}pl 9 25) d ata comparison also showed varied agreement, 
due in part to aerosol inhomogeneity. Again in the 10' 10 to 
10' n m' 1 sr' 1 region. fypsSP, 9.th ^ ala vvere noticeably lower 
due to the undercounting of particles in the 0.8- to 3-/im - 
diameter range. 

In the two transit flight case studies presented here, LOPC 
microphysics supplemented by fig.j f/ftio.59 data indicated that 
a mixed phase of sulfuric acid and ammonium sulfate (acid 
sulfate) for the larger particles was the most likely aerosol 
composition present for most parts of the two flights. These 
particles contribute the most to the infrared aerosol 
backscatter. LOPC data also indicated that smaller particles in 
a subtropical environment are more likely to be externally 
mixed sulfuric acid and ammonium sulfate particles. The GSFC 
lidar data from the very early portion of flight 14 suggested 
the presence of very fine dust and/or sootlike particles in a 
possible pollution plume east of Tokyo, Japan. 

The large continuous set of size distribution data obtained 
with the FSSP sensor enable statistics of infrared backscatter 
to be calculated in the midtropospheric atmosphere over the 
Pacific Ocean. Both the northern and the southern hemisphere 
data sets showed a major backscatter mode occurring around 
3 x 10' 11 m' 1 sr' 1 after a correction is made for large-particle 
underestimation. The mode is in close agreement to that 
reported in earlier work on the presence of an aerosol 
background mode in this region of the atmosphere. 
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